1 



Pseudospin and Deformation-induced Gauge Field in Graphene 

Ken-ichi SASAKI,*) and Riichiro SAITO 

Department of Physics, Tohoku University, Sendai 980-8578, Japan 

The basic properties of 7r-electrons near the Fermi level in graphene are reviewed from a 
point of view of the pseudospin and a gauge field coupling to the pseudospin. The applications 
of the gauge field to the electron-phonon interaction and to the edge states are reported. 

§1. Introduction 

The electronic properties of a single layer of graphite, graphene, x )~ 3 ) have at- 
tracted much attention due to the "relativistic" character of 7r-electrons near the 
Fermi level. The energy band structure of graphene exhibits a linear energy disper- 
sion relation around the two inequivalent, hexagonal corners of the first Brillouin 
zone in the /c-space (the K and K' points). 4 -*' 5 ) The wavefunction (Hamiltonian) of 
7r-electrons has two component (a 2 x 2 matrix) form due to the fact that the unit 
cell of graphene consists of two carbon atoms (A and B atoms). The effective-mass 
Hamiltonian of 7r-electrons around the K point or the K' point is given by linear 
momentum operator, which is relevant to the linear energy dispersion relation of 
graphene. The effective-mass equation is similar to the massless Dirac equation or 
the Weyl equation. 6 ) 

The original Dirac equation for an electron has the form 



d , , , N I —ccr ■ p Imc 



where c is the speed of light, m is the mass of electron, and p (= —iftV) is the 
momentum operator. The Hamiltonian is a 4 x 4 matrix, which is wr itte n in terms 
of the Pauli matrices a = (a x , a y , a z ) and 2x2 identity matrix /. ip(r,t) of 
the Dirac equation is a four-component wavefunction, which naturally explains the 
spin degree of freedom for particle (electron) and the antiparticle (positron) . In the 
massless limit (to = 0), the Dirac equation is split into two equations of a 2 x 2 matrix 
form, that is, the Weyl equations for massless neutrinos. Setting tp = (i^l^r) in 
Eq. (PJ, we obtain the Weyl equation for ipR and ^ as 



d 

ih—^ L (r, t) = -ccr ■ pV>i(r, t), 
d 



(2) 
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We use the Pauli matrices of the form of a x = ( , \ 1 , <?v = ( ? J 1 , and a z — 
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The two component of each Weyl equation represents the spin degree of freedom. 

On the other hand, the effective-mass equations for the K and K' points of 
graphene are written as 

d 

ih—ip K (r,t) = vpo- ■ pip (r,t), 

3 q (3) 
ih—ij K '(r,t) =v F a' -p^ K '(r,t), 

where vp is the Fermi velocity, a = (a x ,a y ), a' = (— a x ,a y ), and p = (p x ,P y )- The 
equation for i/) K of Eq. ([3]) is similar to the Weyl equation for tpR of Eq. ([2]) with 
p z = and by substituting c to vf. The equation for i/j k is similar to the Weyl 
equation for Vl with p z = 0, the substitution of c to vf, and the negative sign 
in front of a y . Although the character of the two component of the effective-mass 
equations of graphene (A and B atoms) and of the Weyl equation (up spin and down 
spin) is different from each other, the equation and the resulting solution for given 
u x and u y are the same. Thus, it is appropriate that the two component structures 
V'K = *C0K' V'k) an d V'K' = t {4 , K'^K') °f the effective-mass equations for graphene 
near the Fermi level is referred to as the pseudospin. A pseudospin structure gives 
a rich variety of interesting physical phenomena of graphene and nanotubes. It is 
known that absence of a backward scattering of an electron in graphene and carbon 
nanotubes is relevant to the nature of pseudospin, 7 - 1 ' 8 ) in which a 2ir rotation of a 
pseudospin wavefunction around the K point in the two-dimensional Brillouin zone 
does not gives the original wavefunction but gives minus sign to the wavefunction. 
The pseudospin in the fc-space behaves similarly to the real spin in the real space. 

The interaction between an electron and an electromagnetic field is given by 
replacing p in the Dirac equation with the kinematical momentum ~k = p — eA(r) 
where — e is the charge of electron and A(r) is a vector potential. The spin of an 
electron is polarized by a magnetic field, B(r) = V x A(r), which can be shown 
explicitly by the Dirac equation. 6 ) In the non-relativistic limit, the Dirac equation 
reduces to the Pauli equation in which the leading interaction between spin and a 
magnetic field is reduced to the Zeeman term: 

-^■B(r). (4) 

Since the equation for graphene is similar to the Dirac equation, the following ques- 
tion arises; What is the field that polarizes the pseudospin? A magnetic field is a 
candidate. The same procedure as that in Eq. @ shows, however, that the pseu- 
dospin is not polarized by a magnetic field because the Zeeman term appears with 
the opposite sign at the K point and at the K' point as 

-^-o--B(r) (K point), 

eh (5) 
+ iLcr.B(r) (K' point). 
2m 

Thus the direction of the pseudospin polarization which is induced at the K point 
by a magnetic field is opposite to that at the K' point by the same magnetic field. 
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Thus the pseudospin of graphene is not polarized by B(r). EE Mathematically, this 
observation leads us to assume a new gauge field A q (r) which has the opposite sign 
of A q (r) at the K point to the K' point as 

d 

ih—ip K (r, t) = v F a • (p + A q (r)) ip K (r, t), 

9 q (8) 
i^ipK'(r,t) =v F a f ■ (p- A q (r))V> K 'M)- 

Then the corresponding field defined by B q (r) = V x A q (r) can polarize the pseu- 
dospin because the Zeeman term appears as the same sign for the K and K' points 
in this case. 

There is an example of the pseudospin-polarized state, that is, a localized state 
appearing near the zigzag edge of graphene, which is called the edge states. 9 ) The 
wavefunction of the edge state has a value only on A atoms when the zigzag edge 
atoms consist of A atoms. Thus we expect that a field B q (r) appears around the edge. 
We will show that the zigzag edge structure is relevant to a field B q (r) which polarizes 
the pseudospin near the zigzag edge. 10 - 1 A gauge field A q (r) and a field B q (r) is 
important not only for the localized edge states but also for the extended states. For 
example, the electron-phonon (el-ph) interaction in graphene can be expressed by 
A q (r), which explains the chirality dependence of the el-ph interaction. n ) In this 
paper we review our studies performed on the gauge field A q (r) by showing that the 
el-ph interaction and the edge boundary are expressed by A q (r). 

Here, we would like to mention the relationship between our work and previously 
published literature on deformation-induced gauge field. Kane and Mele introduced 
a homogeneous deformation-induced gauge field in the effective-mass theory. 12 ) The 
gauge field represents uniform lattice deformations such as uniform bend, twist and 
curvature of a carbon nanotube. The uniform gauge field for the curvature of a 
nanotube changes the boundary condition around a tube axis which is given by a 
generalized Aharanov-Bohm (AB) effect and induces a small energy gap in a chiral 
metallic carbon nanotube. The curvature-induced mini gap was observed by scanning 
tunneling spectroscopy (STS) experiment by Ouyang et al. 13 ' A generalization of the 
gauge field to a local field is necessary for describing a local lattice deformation. In 
the previous paper, 14 ) we generalized the gauge field introduced by Kane and Mele, to 



*' It is noted that even in the massless limit of Eq. @, we obtain a similar coupling term by 
taking the square of the Hamiltonian in the Weyl equation as 

(cct ■ tt) 2 = cl (p - eA(r)) 2 - ehca ■ B(r). (6) 

The first term in the right-hand side of Eq. (|6]) does not concern the spin, but the second term in 
the right-hand side of Eq. (]5J is similar to the Zeeman term which shows that the spin is polarized 
by a magnetic field. 

(vfo- ■ 7r) 2 = Vpl (p — eA(r)) 2 — ehvpa z ■ B z (r), (K point) 

2 * (*0 

{y-pcr' ■ #) = Dpi (p — eA(r)) 2 + ehv\a z ■ B z (r). (K point) 

The opposite sign in front of ehvpcr z ■ B z (r) at the K point and at the K' point shows that the 
pseudospin is not polarized by the magnetic field. 
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include a local lattice deformation of graphene. Further, by defining a deformation- 
induced magnetic field, we explained the local modulation of the energy band gap, 14 ) 
which was observed in the STS measurement for pea-pod (CgO encapsulated carbon 
nanotube) by Lee et al. 15 ^ 

This paper is organized as follows. In §[2l we derive the effective-mass equations 
for graphene and define the pseudospin. We show that the pseudospin determines 
the elastic scattering amplitude of an electron by impurity potentials. In § EJ we 
derive the effective-mass equation in the presence of the lattice deformation, in which 
gauge symmetry for A q (r) is essential to the time-reversal properties of the el-ph 
interaction. In §21 the formulation explained in §[3] is applied to the el-ph interaction 
for optical and acoustic phonon modes. In § [5j we will discuss the edge state by A q (r). 
In § El summary of this paper is given. 

§2. Effective-mass Theory 

In this section, we derive an effective-mass Hamiltonian from the nearest-neighbor 
tight-binding model for a graphene. The nearest-neighbor tight-binding model is 
given by 

^0 = -7oE E + (cf )^+a) , (9) 

ieA a=l,2,3 

where 70 (~ 2.7 eV) is the nearest-neighbor hopping integral, ((cfy) is the 
annihilation (creation) operator of 7r-electron for A-atom at the position rj, and cf +a 
{{cf +a )^) is that for B-atom at r.; +fl (= r.j + R a ) where R a (a = 1,2,3) are vectors 
pointing to the three nearest-neighbor B-atoms from an A-atom (see Fig. [1]). 

We use the Bloch theorem to diagonalize Eq. Q. The Bloch wavefunction with 
wavevector k is defined by 

|^ k > = * £ e *-*(cf)t|o> (, = A,B), (10) 

* u ies 

where the sum on i is taken over the crystal, N u is the number of the hexagonal unit 
cells, and |0) denotes the state of carbon atoms without 7r-electrons. The off-diagonal 
matrix element of TCq is given by 

<3$|Wo|s£> = -70 E /a(k) = - 7 o/(k), 

(11) 

<#bWa> = -7b E /a(k)* = -7o/(k)*, 

a=l,2,3 

where / a (k) = e lk Ra and /(k) is given by 16 ' 



/(k)= /«( k )- 

a=l,2,3 



(12) 
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Fig. 1. A hexagonal unit cell of graphene consists of A (solid circle) and B (open circle) atoms, ai 
and &2 are the unit vectors (a denotes the lattice constant: a = |a;| = V3a cc where a cc is the C- 
C bond length.). R a (a = 1, 2, 3) are vectors pointing to the nearest-neighbor B atoms from an 
A atom. For the xyz coordinate system, R a are written as Ri = a cc e y , R2 = — (\/3/2)a cc e x — 
(l/2)a cc e y , and R3 = (y/3/2)a cc e x — (l/2)o C ce B where (e y ) is the dimensionless unit vector for 
the a;-axis (j/-axis). Local modulations of the hopping integral are defined by S'yg (ri) (a = 1, 2, 3) 
where ri (i £ A) is the position of an A-atom. 



The diagonal matrix element of 7io, (^(Wol^), (s =A,B) can be taken to be zero. 
The energy eigenequation is written by a 2 x 2 matrix form as 

*U -.„( ° W 4 '-' 



The energy band structure of graphene is obtained by solving 

det ( ^ (k) 7o/(k) "l - (14) 

whose solution is E(k) = +7o|/(k)|, (— 7o|/(k)|) for the conduction (valence) energy 
band. The conduction energy band and the valence energy band touch with each 
other at the K point, k F (= (4-7r/3a, 0)), and at the K' point, — k F , where |/(k)| 
vanishes. By expanding / a (k) in Eq. (jlip around the wavevector of k F (the K 
point), we obtain 

/ a (k F + k) = / a (k F ) + i/ a (k F )k • R a + • • • • (15) 



Here k (= (k x ,k y )) is measured from the K point. Using k F = (47r/3a, 0), we get 
/ x (k F ) = 1, / 2 (k F ) = e~ i2 -i , and / 3 (k F ) = e +i ^ . Substituting / (k F ) into Eq. (p|), 
we obtain from Eq. (Ilip that 

K F+k |Wo|^ B kp+k ) = 7o^(fc, - iky) + 0(k% 

3a ^ 
(^B F+k |^o|^A F+k ) = 7o^(^ + iky) + 0(k 2 ), 
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where we used (^ F \H \^ F ) = -7o/(k F ) = 0. 
From Eq. (fTUj) . we see that Eq. (fTBj) becomes 



E(k F + k 




37 a c 








(17) 



By introducing the Fermi velocity as vp = 7o3a cc /27i (~ 10 6 m/s), the momentum 
operator p = —ihX7, and the Pauli matrix a = (a x , a y ), we obtain the effective-mass 
Hamiltonian near the K point as 

/K 



H 



o 



Vp(T ■ p. 



(18) 



7^0 is a 2 x 2 matrix which operates on the two component wavefunction: ^ K (r) = 
*(V*A( r )> V*b 00) where V ; A( r ) anc ^ V'b( 1 ') are ^ ne Bloch wavefunction for 7r-electrons 
of the A and B atom in the unit cell. By introducing 0(k) which is defined 
by an angle of k = (k x ,k y ) measured from the k x -axis, then we 
|k|(cos0(k),sin@(k)), and 



Ti 





-i<9(k) 



-iO(k) 




(19) 



The energy eigenvalue of Eq. (|19p is given by ±t>p|p| and the energy dispersion 
relation shows a linear dispersion relation near the K point as is known as the Dirac 
cone. At the K (or K') point, the valence and conduction bands touch to each 
other and thus we call the degenerated point as the Dirac point. The eigenstates for 
E = +vp\p\ and E = — i>f|p| are given by 



V> c K k (r) 



-i©(k)/2 
-i©(k)/2 



V£k(r) 



ikr 



e -i0(k)/2 \ 

-e 



(20) 



/25 W/ V ' V25 

which are a conduction band state and a valence band state with the wavevector 
In Eq. ([2"0"j) . 5 denotes the surface area of graphene. In the case of single wall 



**) 



We note that 



^c K k (r) 



1 

,+;e(k) 



+»e(k) I i 



(21) 



are also eigenstates of Eq. (|18[) with an additional phase factor of e ±e ' k ' //2 . These wavefunctions 
do not change their values under a 2n rotation in the fe-space, <9(k) — > <9(k) + 27r. This behavior is 
different from Eq. (|20|l , in which the wavefunction of Eq. (|20[l changes the sign after the 27r rotation. 
However, the Berry's phase for Eq. (|21|l which is defined by 



d0 



1 

71 



l,e 



-ie(k) 



1 

71 



1 

,+ie(k) 



(22) 



gives an extra phase shift of it. On the other hand, the Berry's phase for Eq. (I20|l vanishes because 



(10 



1 

71 



e +ie(k)/2 )6 - 



te(k)/2\ 







/-ie(k)/2y 






e +;e(k)/2 J 



0. 



(23) 



Thus, a 2-7T rotation of a pseudospin wavefunction around the K point in the two-dimensional Bril- 
louin zone gives minus sign to the wavefunction in the both cases of Eqs. (|20|l and (|2ip . It is 
convenient to use the wavefunction of Eq. ((20} , since the effect of Berry's phase is included in the 
wavefunction. 
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carbon nanotube (SWNT), the k x -axis is taken in the circumferential direction on 
the cylindrical surface and the k y -axis is defined by the direction of a zigzag nanotube 
axis (see the coordinate system in Fig. QJa)). We see in Eq. (|20p that the energy 
eigenstate for the valence band, ^f k ( r ) * s g iven by a z ip^ k (r). This is because of a 
particle-hole symmetry of the Hamiltonian: TCq(j z = —OvI-Oq. PI 

Similarly, the effective-mass Hamiltonian for the K' point is given by expanding 
/o(k) around — kp (the K' point) in Eq. (fTT|) as 

Hf = v F a' • p, (24) 

where a' = (—cr x , a y ) and 7{q' operates on a two-component wavefunction: tp K ' (r) = 
*("0a ( r )' V'b ( r ))- The energy eigenstates for the conduction and valence energy band 
are given, respectively, by 

e *'-r / e +i0(k')/2 \ gik'-r / +i0(k')/2\ 

^ck'(r) = ^= ( e -ie (k ')/2 J . ^,k'(r) = a^ c ,k'(r) = -j= L-ie(k')/2 J ( 25 ) 

where k' — ikl, = IkHe - * ^'). 

The linear energy dispersion relations near the K and K' points, E = ±vp|p|, are 
contrasted to the non-relativistic energy dispersion relation of |p| 2 /2m where m is 
the effective-mass of the particle. The effective-mass for graphene can be understood 
to be zero from the definition of relativistic energy E = iy^c 2 + m 2 c 4 = ±c|p| 
for m = 0, where c is substituted to i>f (~ c/300). The wavefunction with two 
components is defined by the "pseudospin" . The pseudospin up (down) state *(1,0) 
(*(0, 1)) corresponds to the wavefunction which has a value only on A (B) atoms. As 
we see in Eq. (|20p . if we rotate k by 2tt around the K point k = 0, that is, (9(k) + 2ir, 
then V ; J C @(k)+27r( r ) = — ^c C ©(k)( r )' wn ich is the same structure for a real spin under 
a 2tt rotation in the real space. 

Here, we give an example to show that the pseudospin is relevant to a vanishing 
matrix element for the lowest order backscattering amplitude. 17 ** When a carrier in 
the conduction band is denoted by | kc) , then the matrix element of the backscattering 
process is written as (— kc|Wi mp |kc). When the impurity potential is long-range 
compared with the lattice constant, the potential is modeled by a diagonal form and 
the matrix element becomes 

{ kc\n imp \kc) - j 2S e y e -ie(- k y 2 J y Q y( r )^ e +«e(k)/2^ - u ' W 

where we use Eq. (|20p and Q{— k) = (9(k) ± tt. This means that the interference 
between the two component pseudospin makes the back scattering matrix element 
vanish. In general, there are many impurities in a sample, which gives rise to the An- 
derson localization in the disordered systems. There is an interesting way to explain 
the ballistic transport of carbon nanotube even in the many scattering events. ^ For 

*' Ti^a z = —a z ?ij§ is obtained directly by Eqs. (|18[) and a z Ci = —OiO z where i — x or y. Then 
K(a z ^ k (r)) = -ff,?#< k (r) = -a^ F |p|^ K k (r) = (-v F \p\)a z ^ k (r). 
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a back scattered wave, we can find a time-reversal scattering wave whose wavefunc- 
tion has an additional phase shift of vr. The time reversal pair of scattered waves 
cancel with each other, which leads to the absence of the backward scattering. 7 ) 
This is contrasted with the standard concept that impurity scattering gives rise to 
the Anderson localization for low-dimensional systems. 



§3. Gauge Field for a Deformed Graphene 



The dynamics of the conducting electrons in graphene materials are different 
from those of ideal flat graphene, because in the former case, there are shape fluc- 
tuations (cylindrical shape, phonon vibration, etc.) that result in the modification 
of the overlap matrix elements, of nearest-neighbor 7r-orbitals and of the on-site 
potential energy. We refer the modification of the nearest-neighbor hopping integral 
as the off-site interaction and a shift of the on-site potential energy as the on-site 
interaction. The modification of the effective-mass equations by the off-site (on-site) 
interaction is discussed in Sec. 13.11 (Sec. I3.2|) . We discuss time-reversal symmetry 
property of the effective-mass equations including the off-site and on-site interac- 
tions in Sec. 13.31 

3.1. Off- site interaction 

First we consider the perturbation from the off-site interaction in which only 
off-diagonal matrix element has a non-zero value. A lattice deformation induces a 
local modification of the nearest-neighbor hopping integral as —70 — > —70 + #7g(ri) 
(a = 1, 2, 3) (see Fig. [[J. We define the perturbation TCi as 

Wi = £ E WO (( c ?+a) f cf + (cf ) f cf+a) • (27) 

iSAa=l,2,3 

The off-site matrix element of Ti\ with respect to the Bloch wave functions in Eq. (|10p 
with k and k + <5k is given by 

(SPf^lWilffg) = ±- E E 57o a (r,)/a(k)e*, 

% e Aa=l,2,3 (28) 

<^ +5k |Wi|^> = — E E 57o(ri)/a(k)*e- 5k -^ +R "). 

u i G Aa=l,2,3 

Here we consider the two possible cases for 5k. When <5k is sufficiently small com- 
pared with the reciprocal lattice vector, k near the K (or K') point is scattered to the 
k' = k + 5k within the region near the K (or K') point, which we call the intravalley 
scattering. On the other hand, if 5k is comparable to the distance between the K 
and K' points, we expect scattering from K (K') point to K' (K) point, which we 
call the intervalley scattering. 



*' Here the time reversal pair of scattering is defined for the scatterings within one valley of 
the K (K') point. Later, we will use "time reversal symmetry" for the pair of k around the K point 
and k' around the K' point, when we consider the intervalley scattering. 
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3.1.1. Intravalley Scattering 

In the case of intravalley scattering, when k is measured from k F , we get 

(^^I^I^V^E E ha(r t )f a (k F )e- iS ^ + 0(5k5 la ), 

u ieA a=l,2,3 

(a/ k F+ k + <5k |Wi| ^k P+ k ) = 1 ( 5 7a (r i )/a(k F )*e-' 5k -^ + 0(Skd la ), 

u ieA a=l,2,3 



(29) 



The correction indicated by 0(5k5j a ) in Eq. f|29|) is negligible when |<5k| <C |kf 
Substituting /i(k F ) = 1, / 2 (k F ) = e~ i2 r and / 3 (k F ) = e +i2 r into Eq.(Tj9D, we get 



(30) 



^k F+M|Wi| ^k F+k) = vf { ^ (r . } _ a q (r . )} e - l5k . r% 

where A q (r) = (A%(r),A%(r)) is defined by ^(r) (a = 1,2,3) as 
UF ^(r) = *rf(r) - ~ (5 7o 2 (r) + <5 7o 3 (r)) , 

^(r) = ^(5 7 o 2 (r)-5 7 o 3 (r)). 
Since the diagonal term vanishes, that is, (&*\Hi\&¥) = (s = A,B), Eq. (J30 



(31) 



shows that 7i\ is expressed by t> F <x • A q (r) in the effective-mass Hamiltonian. Thus, 
the total Hamiltonian of a deformed graphene is expressed by 

nf+Hf = ^ F <x-( P + A q (r)). (32) 

Similarly, by expanding / a (k) in Eq. ([28]) around — k F of the K' point, we find Ji\ 
appears as — w f <t' • A q (r). Thus, we obtain the effective-mass Hamiltonian for the 
K' point as 

nf' + nf = v F <r'-(p- A q (r)). (33) 

The corresponding Schrodinger equation for Eqs f|32|) and (|33p are given in Eq. (J3]). 14 ) 
Thus, the off-site interaction can be included in the effective-mass equations as 
a gauge field, A q (r). We call A q (r) the deformation-induced gauge field and 
distinguish it from the electromagnetic gauge field A(r). 14 ) TLq + 7i\ does not break 
time-reversal symmetry even tough the A q (r) appears as a gauge field, because the 
sign in front of A q (r) is opposite to each other for the K and K' points. This contrasts 
with the fact that A(r) (magnetic field) violates time-reversal symmetry because the 
sign in front of A(r) is the same for the K and K' points as p — > p — eA(r). The 
time-reversal symmetry will be discussed in detail in Sec. 13.31 The rotation of A q (r) 
or the deformation-induced magnetic field is defined by B q (r) = V x A(r) as 

B?(r) ^_^. 
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The direction of B z (r) is perpendicular to the graphene plane B q (r) = (0, 0, B% (r)). 
The deformation-induced magnetic field changes the energy band structure, which 
will be shown in § [5j The B z (r) couples with the pseudospin through a z term. This 
is shown by taking the square of the effective-mass equations, as we see in Eqs. ([6]) 
and ©, 

E^ K (r) =v 2 F {(p + A^(r)) 2 + hBj(r)a z } V> K (r), 

E^ K '(r) = vl {(p - A^r)) 2 + hB*(r)a z } ^ K '(r), ( '" ) 

where (r)cr 2 terms modify the direction of the pseudospin to the direction opposite 
to i?2 (r) in order to decrease the energy for the both K and K' points. Here let us 
stress a notable difference between B z (r) and the usual magnetic field, B z (r). For 
Bz(r), a uniform magnetic field along the z-axis over a sample is possible to realize 
the uniform spin polarization (the Zeeman term). However, since B z (r) is generated 
by the modification of 70, the field B z (r) appears only near the defect region with- 
out periodic symmetry such as the zigzag edge 10 ) which gives local polarization of 
pseudospin. In general, when a graphene sample contains an equal number of A and 
B atoms to each other, we have 

I B^(r)d 2 r = 0, (36) 
Js 

where S denotes an integral for a whole sample. 

3.1.2. Intervalley Scattering 

Next, we consider the matrix element of 7i\ for the intervalley scattering. Let us 
define the wave vector of a Bloch state near the K point; kp + k and that of another 
Bloch state near the K' point; — kp + k' (|k|, |k'| <C | kjr | ) , then the matrix element 
from the K' to K points is written by using Eq. ()28[) as 



u ieAa=l,2,3 

^ArE E W0/a(kF)*e- i2kF - r V( k '- k > r < 

u ieA a=l,2,3 

= W E + iA v^} e" l2k - r V( k '- k )- r % (37) 

V A 



and 



(^ kF+k |Hi|^ kF+k ') = ^E E 57o(r,)/a(k F )*e- l2kF -^e l ( k '- k )-^e- 4k - R « 

u ieA 0=1,2,3 

^o a (ri)/a(k F )*e- i2k -^e^ k '- k )-^ 

u ieA 0=1,2,3 

= W E K( r *) + M "( r ^} e" i2k - r V( k '- k K (38) 
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In Eqs. (J37J) and (138D . from the first line of the right-hand side to the second line, 
we approximate e lk R<1 ~ 1 and e~* kRa 1 since we consider the case of |k'| <C |kp| 
and |k| <C |kp|. Equations (|37p and (|38|) show that the deformation-induced gauge 
field A q (r) (Eq. (|3ip ) can be applied to the matrix element for intervalley scattering, 
too. It is noted that Eqs. ([37]) and (f38j) are identical to each other up to the leading 
term. The matrix elements from the K point to the K' point are given by replacing 
k F by -k F in Eqs. ([33) and ([38]) as 



(39) 



^-k F+ k' |Wi| ^k P +k ) = j2 | A q( ri ) _ iA^Ti)} e +i2k F T ie -i(k'-k). rj 



ieA 



ieA 



In the effective-mass model, the intervalley scattering appears as the off-diagonal 
terms in a 4 x 4 matrix: 

/^ K (r)\ /<r-(p + A q (r)) ^(r)(7« \ /^ K (r)\ ^ 

U K '(r)J " V <r' • (p - A%r))) {^(r)) ' (40) 

where we define a deformation-induced complex scalar field 4>(r) 

^(r)^(^(r) + ^(r))e- 2k - r . (41) 

When A q (r) is a smooth function of r, the effect of intervalley scattering by the 
complex scalar field </> q (r) can be neglected because q (r) becomes a rapidly oscillat- 
ing function of r due to the factor e~* 2kF ' r , and the value of matrix element between 
the K and K' points becomes negligible, too. On the other hand, when A q (r) is 
a rapidly oscillating function of r, </> q (r) becomes a smooth function of r and the 
matrix element between the K and K' points is not negligible. For example, phonon 
modes at the K point contribute to the intervalley scattering. Since the vibrational 
amplitude of the K point phonon modes are proportional to e~* kp ' r ', A q (r,) has the 
same periodicity of e +2ikF ' ri where we use the fact e +3lkF ' ri = 1. Then, <^> q (rj) would 
behave as a constant in the real space. 

3.2. On-site interaction 

Now we consider the on-site interaction by a defect of the crystal. A lattice 
deformation gives rise to not only a change of the transfer integral between A and B 
atoms but also a change of the potential at the A (B) atom <pA (^b) which we call 
the off-site and on-site deformation potential, respectively. We denote the on-site 
deformation potential by a 2 x 2 matrix for tp = t (tp K , ip K ) 

Let the displacement vector of A-atom at r-j is ua(i*j) and that of B-atom at rj is 
ue(rj), then the deformation potential of A-atom at r^, (r-j), is induced by the 
relative displacements of three nearest neighbor B-atoms from the A-atom (113(1*, + 
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R ) - u A (rj)) as 



A (r;) = j^- ^2 R a • (u B (rj + Ro) - u A (r;)) , 

acc 0=1,2,3 



(43) 



where g on denotes gradient of the atomic potential at rj, and t denotes 3a cc /2. Here 
we assume that |u B (rj + R a ) — u A (rj)| "C a cc and that 4>a(^i) depends linearly on 
the relative displacement vector. 




B uif(ri^|R 3 ) 
A 



Fig. 2. Displacements of B-atoms at Vi + R a (a = 1, 2, 3), ub (r^ + Ra), give rise to a deformation 
potential at A-atom of ri. 

By expanding u B (ri + R 2 ) as u B (r-j + R 2 ) = u B (rj + Ri) + ((R 2 - Ri) • V)u B (r; + 

Ri) H and u B (rj + R 3 ) as u B (rj + R 3 ) = u B (rj + R x ) + ((R 3 - Ri) • V)u B (r; + 

Ri) + • • • , we see that Eq. ([4*3]) can be approximated by 

#A(r<) =5onV-u B (r, ; + R 1 ) + --- , (44) 

where we have used Yla=i 2 3 = 0- It is noted that a general expression for the 
deformation potential, Eq. (|44p . is valid in the case that u B (r) is a smooth function 
of r. When this is not the case, we have to use Eq. (|43p . In the continuous limit, we 
may use r to represent the positions of both A and B atoms in the unit cell, then 
we have 4>a ( r ) = g nV • u B (r) + • • • . Similarly, the deformation potential of B-site 
of ^ + Ri is given by 

B (r 4 + R 1 ) = |^ -Ra-(u A (r l + Ri-R a )-u B (r i + R 1 )). (45) 

tClcc a=l,2,3 

By using u A (rj + Ri - R 2 ) = u A (r^ + ((Ri - R 2 ) • V)u A (rj) H and u A (ri + 

Ri — R3) = u A (rj) + ((Ri — R3) • V)u A (rj) + • • • , we see that Eq. (|45l) can be 
approximated by 

B ( ri + Ri) = 5on V-u A (ri) + --- . (46) 

Thus, for the intravalley scattering, we may rewrite Eq. (f4"2|) using Eqs. ([4*4"]) and 
(USD as 



( 'V-u B (r) 
n ° n ~ 9on 1 V • u A (r) ' " 



(47) 
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According to the result of density- functional theory by Porezag et a/., 18 )' 19 ) we will 
use the parameter for g on (=17eV). For the later discussion of el-ph interaction of 
acoustic 

s(r) S "aM + ubM (48) 

and optical 

u(r) = u B (r)-u A (r), (49) 
phonon modes, we rewrite Eq. (|47p using the Pauli matrices as 

n on = 9 -flV • (u A (r) + u B (r)) + 9 -fa z V ■ (u B (r) - u A (r)). (50) 

Comparing with the expressions for A q (r) and ^ q (r) in Eqs. (|31|) and (|41|) . and for 
Eq. (|50p . the effective-mass Hamiltonian of the defect of the crystal is given by 

\VF0 q (r)V x W K ' J' (51) 

where we define 

H K = v F * • (p + A q (r)) + 5on /V • s(r) + ^a z V • u(r), 

(52) 

H K ' = v F a' ■ (p - A q (r)) + g on IV ■ s(r) + ^a z V • u(r). 

3.3. Time-reversal and pseudospin symmetries 

Here we show that A q (r) keeps time-reversal symmetry of graphene system. The 
time-reversal operation, T, is defined by exchanging the K point and the K' point, 
and taking a complex conjugation of the four component wavefunction as 



and K is the complex conjugate operator. In order to check that Eq. (|53|) is time- 
reversal operation, it is useful to introduce the vector potential for electro-magnetism, 
A(r), into the effective-mass Hamiltonian as 

/<x-(p + A q (r)-eA(r)) ^(r)a x \ 

n ~ F V <A q (r) V, a» ■ (p - A q (r) - eA(r)) J ' ^ 

Then the electromagnetic current operator, J(r), is given by 

J(r) = = ~ eVF (o £)" (55) 

Using T in Eq. (|53p . we can show that 

T- 1 J(r)T = -J(r), (56) 
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where we have used K~ 1 crK = —cr' and K^er'K = —cr. The negative sign on the 
right-hand side of Eq. ([56]) shows that Eq. ([53]) is time-reversal operation. 
When we apply T to the wavefunction of Eq. (|40p , we get 



i> K \r)*\ _ Ar-(P + Ai(r)) <j>%r)a x \ 



r 



S U K ( r )V~ VF V WV, ' < r'-(p-A ( i(r))J ^ K (r)*J- (5?) 
By taking the complex conjugation (c.c.) of the first row of Eq. (|57p . we get 

E^ K '{r)* = v F cr • (p + A q (r))^ K '(r)* + ^(r^^r)* 
^ ^ K '( r ) = upcr' • (p - A q (r))^ K '(r) + v F (^{r)* a x ^ K {v) (58) 

where we have used cr* = —cr' and p* = —p. Similarly, by taking the complex 
conjugation of the second row of Eq. (|57p . we get 



i# K (r)* = v F 4>\r)*a x ^ K \v)* + v F cr' • (p - A q (r))^ K (r)* 

Eip K (r) = v F ^{v)a x ^ K ' {r) + v F er • (p + A q (r))^ K (r). (59) 



The equations in the last lines of Eqs. (|58p and (|59p are nothing but the second and 
first row of Eq. (|40p . respectively. Therefore, Eq. (|54p (or Eq. (|40p ) is symmetric 
under the time-reversal transformation (when A(r) = 0). 

When t/> q (r) = 0, it is useful to define an operation S that transforms the electron 
of k to that of — k within the same valley. The operation S is defined by 

^ K (r)WM K (rn ^_c_^ 



^ (r) J - \a y ^ (r)*J ' whare 5 = (,0 /J ^ (60) 

This operation is referred to as effective time-reversal symmetry or special time- 
reversal symmetry because we obtain S~ 1 3(r)S = — J(r) for the special case of 
</> q (r) = 0, A q (r) = 0, and A(r) = 0. By applying Eq. (|6U|) to the wavefunction of 
Eq. (|40p . we get for the first row as 

Ea y ^ K (v)* = v F cr • (p + A q (r))<7^ K (r)* + v F ^(r)a x a y ip K ' (r)* 
% E^ K {v)* = v F cr' • ( P + A q (r))^ K (r)* - v F ^(r)a x ^ K ' (r)* 
C A E^ K (r) = v F cr • (p - A q (r))^ K (r) - v F ^(r)*a x ip K '(r) (61) 

Here, from the first line to the second line, we multiplied a y to both sides and used 
Oycrcjy = cr' . To get the last line, we took the complex conjugation of the second line 
and used (cr')* = —cr and p* = —p. The special operation S becomes a symmetry 
when A q (r) = 0. Note that A q (r) = results in c/> q (r) = because of Eq. (JIT]). 

We comment on the definition of the time-reversal operation and the on-site 
interaction. The transformation defined by Eq. (|53|) does not contain any Pauli spin 
matrix since time-reversal symmetry has nothing to do with the pseudospin degree 
of freedom. It is also noted that the pseudospin degree of freedom (a y in Eq. ([60]) ) 
is necessary to define the special symmetry of Eq. (|60p . Even in the presence of the 
on-site deformation potential of Eq. (|42p , the time-reversal symmetry of Eq. (|53p is 
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valid since Eq. (|53j) does not contain the Pauli spin matrix. However, the special 
symmetry of Eq. (|60p is lost when Eq. (|42p is not symmetric about the sublattice. 20 ^ 
Namely, when we take into account a non- vanishing a z term of Eq. (|50p for Eq. (|6ip , 
the symmetry of Eq. (|60p is broken even when A q (r) = 0. 

§4. Electron-Phonon Interaction 

In this section, we apply the A q (r) for phonon mode to describe the el-ph interac- 
tion. There are six phonon modes in graphene: in-plane longitudinal optical/acoustic 
mode (LO/LA), in-plane tangential optical/acoustic mode (iTO/iTA), and out-of- 
plane tangential optical/acoustic mode (oTO/oTA). 16 ) We consider long wavelength 
in-plane optical phonon modes; the LO and iTO phonon modes near the r point. 
Using Eqs. (|3ip and (|50p . we will derive A q (r) and 7i on for the el-ph interaction 
of the LO/iTO modes. The material in this section has been used to analyze the 
LO/iTO phonon frequency shift in metallic SWNTs. 11 ) The phonon frequency shift 
is observed as a function of the Fermi energy (the Kohn anomaly) in the Raman 
spectroscopy experiments. 21 )~ 26 ) 

4.1. r point optical phonon modes 

The el-ph interaction originates from a change of the atomic potential due to 
a vibration of a carbon atom. In Fig. [3l we show a change of the atomic potential 
whereby <57o(r) is induced. The atomic potential whose origin is located at r is 
denoted by solid curves and the shifted potential whose origin is located at r + Ar is 
plotted by dashed curves. The deformation potential is given by <P(r + Ar) — ^(r) = 
Ar ■ V<P(r) + • • • , which gives rise to matrix element between the 7r-electron at A- 
atom and the 7r-electron at B-atom as (^(r + R a )|V<P(r)|<pA.(r)) • (R a /a cc ). In this 
paper we denote this matrix element as g g/£, that is, (^(r + R-a)|V#(r)|<P r A(i")) • 
(R a /a cc ) = g Q ff/£- According to density functional calculation by Porezag et al., 18 ^ 
we have g oS /£ = 3.0eV/A 19 )' 27 ) where £ = 3a cc /2 = 2.13A. 




Fig. 3. The atomic potential of a carbon atom at r is denoted by solid curves and the shifted 
potential of the carbon atom at r + Ar is plotted by dashed curves. The deformation potential 
at the B atom is given by <P(r+Ar)—$(r) — Ar- V<£(r)+- • • . The deformation potential gives rise 
to a scattering of a 7r-electron from A-atom to the B-atom as (^B(r-|-R a )|V<?(r)|S'A(r))-(R.a/ £ icc)- 

When Ar = u(r) (u(r) is defined in Eq. (|49p ). the change of the C-C bond- 
length is given by u(r) • (R a /a cc ). Thus 57g(r) for the LO and TO modes is given 
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by 



*tf(r) = ^u(r) ■ ?S. (62) 



By putting Eq. (f6"2"j) into Eq. ([3"Tj) . we obtain 



UF^(r) = — u(r) ■ Ri - 



2 

WF ^(r) = gL^u(r).(R 2 -R 3 ) 



(63) 



Because Ri — (R2 + R3V2 = te^ and \/3/2(R2 — R3) = — £e x (see the caption of 
Fig. [I]), we can rewrite Eq. ([63]) by A q (r) as 

v F (A%(r),A«(r)) = M (%(r) _^ (r))) (64 ) 

where u x .(r) = u(r) • e x and %(r) = u(r) • e^. 28 ) For the on-site el-ph interaction of 
Eq. (|50p . we can neglect the term which is proportional to 00 since the LO and TO 
optical phonon modes satisfy ua(i") = — ue(r). Then, we get 

n on = ^a z V-u(r). (65) 

The resulting effective-mass Hamiltonian for the K point and K' point become 

H K = v F a ■ (p + A^(r)) + ^a z V ■ u(r), 

n ( 66 ) 
H K ' = v F a' • (p - A») + ^a z V • u(r), 

where A q (r) is given by Eq. (|64p . Putting Eq. (|64p into Eq. (|34|) . we see in Eq. (|66p 
that V • u(r) is proportional to the deformation induced magnetic field since 

V-u(r) = -^^(r). (67) 
9oS 

Thus, if Bz (r) 7^ 0, the energy spectrum for an electron has a local energy gap around 
E F because the terms proportional to a z in the Hamiltonian Eq. (|66p work as time 
dependent mass term. l!j However, the optical phonon mode does not generate a 
static (time independent) energy gap because u(r) is oscillating as a function of time 
as u(r) cos(wt) so that the time average of B z {r) field vanishes. However, the energy 
gap can be oscillating as a function of time when the vibration of atoms is coherent 
in the space. 

Next we consider an out-of-plane TO (oTO) phonon mode. The oTO phonon 
eigenvector is pointing in the direction of perpendicular to the 2D plane and does 
not give rise to a first order, in-plane bond-length change. Thus, the off-site and 
on-site el-ph coupling for the oTO mode are negligible as compared with those for 
in-plane LO and TO modes. 

The mass term is defined by mo z in the effective-mass Hamiltonian where m produces an 
energy gap in the energy spectrum. 
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4.1.1. Phonon softening for q = LO/iTO phonon mode 

For a uniform (u x ,u y ) in the real space, that is, for the LO and TO modes 
with the phonon wave vector q = 0, the el-ph interaction for the LO and TO 
phonon modes is given only by A q in Eq. (|66p since V • u(r) = 0. For the r 
point optical phonon modes, we consider a virtual electron-hole pair creation by the 
el-ph interaction, which contributes to the phonon frequency shift in second-order 
perturbation theory (see Fig. Sfa)). The phonon frequency shift as a function of 
the Fermi energy has been observed in Raman spectroscopy for graphene 21 )' 22 ) and 
metallic nanotubes. 23 ^ 26 - 1 

The el-ph matrix element for the electron-hole pair generation near the K point 
is given from Eqs. (|20p and (|64p by 

(eh(k)\Hf\u) = f (V c K k (r))^ F <x- A°PV4(r)d 2 r 

9oS (e +t2 ¥ L \ ( u y + iu x \ ( e^ 2 ^ 1 \ . , 

" 2^ { e -tep ) [uy -iu x o ) [_ e+ i^ ) ■ (68) 

Let us introduce the angle a between vector u = (u x ,u y ) and the x-axis, then by 
putting u x = ucosa and u y = usina (u x + iu y = ue ia ) into Eq. ([68]) we obtain 

(eh(k)\Hf\u a ) = -ig oS — cos(0(k) - a). (69) 

^cc 

The el-ph matrix element for the electron-hole pair generation near the K' point, 
(eh(k)|?i^ \uJ a ), is the same as Eq. (j69j) . In Eq. (|69|) . a depends on the LO and TO 
modes in the case of SWNTs, while a can not be defined uniquely in the case of flat 
graphene. Let us consider the case of a zigzag SWNT. Then, we denote x (y) as a 
coordinate around (along) the axis as shown in Fig. QJb). In this case, u = (u x ,0) 
(a = 0) is assigned to the TO phonon mode while u = (0, u y ) (a = vr/2) is assigned 
to the LO phonon mode. By calculating Eq. (|69p for the TO mode with a = and 
for the LO mode with a = vr/2, we get 

(eh(k)\Hf |o"ro) = -ig s — cos@(k), 

(70) 

(eh(k)|ftf \oj lo ) = -ig oS —smO(k). 

For a "metallic" zigzag SWNT without the curvature effect, 29 )' 30 - 1 we obtain 
(9(k) = ±7r/2 for a metallic energy band with k x = (see Fig. 2Kb)). Then, Eq. ([TO]) 
tells us that only the LO mode gives rise to a non-vanishing el-ph matrix element and 
the TO mode does not contribute to the electron-hole pair creation. The amplitude 
for an electron-hole pair creation depends strongly on the curvature effect which 
shifts the relative position of the cutting line from k x = (see Fig. 2Kb)). When 
the curvature effect is taken into account, cos@(k) = k x /(k x + k y ) 1 ^ 2 is nonzero due 
to k x ^ 0. Thus, the low energy electron-hole pair can couple to the TO phonon 
mode and the matrix element becomes maximum when k y = 0. On the other hand, 
the high energy electron-hole pair still decouples to the TO phonon mode since 
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cos(9(k) — * for \k y \ 3> \k x \. This leads to a phonon frequency hardening for the 
TO phonon mode of a zigzag nanotube when the Fermi energy is located near the 
Dirac point. n ) 




Fig. 4. (a) An electron-hole pair on the cutting line near the K point is excited by the r point 
phonon mode (q = 0). The intermediate electron-hole pair state contributes to the phonon 
frequency shift, (b) The k x (k y ) axis corresponds to the zigzag SWNT circumferential (axis) 
direction. The amplitude for an electron-hole pair creation (Eq. (J70j ) depends strongly on the 
relative position of the cutting line from the K point. If the cutting line crosses the K point, 
then the angle <9(k) (= a.rcta.n(k y / k x )) takes n/2 (— tt/2) values for k y > (ky < 0). In this 
case, the LO mode strongly couples to an electron-hole pair, while the TO mode is decoupled 
from the electron-hole pair. 



The description of deformation induced gauge field for a lattice deformation 
(Eq. (|3ip ) is useful to show the appearance of the curvature- induced mini energy 
gap in metallic zigzag SWNTs. 12 ) For a zigzag SWNT, we have = and 5jq = 
5jq from the rotational symmetry around the tube axis (see Fig. [1]). Then, 
using Eqs. (|31~]). we get that Ax = —S^q/vf and Ay = as the curvature effect. The 
cutting line of k x = for the metallic zigzag nanotube is shifted by a finite constant 
value of Ax because of the AB effect for the deformation- induced gauge field A q (r). 
This explains the appearance of the curvature-induced mini energy gap, 

60(meV-nm 2 ) 
-fcgap ^ ' 'y' 1 ) 

where dt is the diameter of a metallic zigzag SWNT. 
4.1.2. 0, LO/iTO mode 

The description of the el-ph interaction as a gauge field A q (r) can be extended 
to show the decoupling between the TO mode with q ^ and the electrons using a 
gauge symmetry argument. n ) The TO phonon mode with q / does not change 
the area of the hexagonal lattice but instead gives rise to a shear deformation for the 
hexagonal unit cell of graphene. Thus, the TO mode (uxo( r )) satisfies 

V-u TO (r) = 0, Vxu TO (r)/0. (72) 

Thus, the on-site deformation potential is zero in Eq. ([66]) . Using Eqs. ([64"}) and (fT2j) . 
we see that the deformation-induced magnetic field, Bg(r), becomes zero instead the 
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divergence of A q (r) exists, which are 



(V x A q (r)) • e 2 = B*(r) 



9oflf 



V • UTo(r) = 0, 



VFa cc 



(73) 



V • A q (r) 



9oS 



(Vxu TO (r))-e z /0. 



VFCl, 



cc 



In this case, A q (r) can be represented by the gradient of a scalar function, $a( r )> 
as A q (r) = V^afr). Since we can choose a A q (r) = gauge in Eq. ([66]) by a 
redefinition of the phase of the wavefunction as ■0 K (r) — > exp(— il^ a (r)//i)^ K (r) and 
V> K '(r) — * exp(+i'P r a (r)//i)V' K '(r), 14 ) the A q (r) field in Eq. ([66]) disappears for the 
TO mode with q / 0. This explains why the TO mode with q / completely 
decouples from the electrons and that only the TO mode with q = couples with 
electrons. In this sense, the TO phonon mode at the r point is anomalous since 
the el-ph interaction for the TO mode can not be eliminated by a phase of the 
wavefunction. On the other hand, the LO phonon mode with q / changes the 
area of the hexagonal unit cell while it does not give rise to a shear deformation. 
Thus, the LO mode (uLo(r)) satisfies 



Using Eqs. (|64p and f|T4j) . we see that the LO mode gives rise to a deformation- 
induced magnetic field as 



Since a magnetic field changes the energy band structure of electrons through the 
mass term, the LO mode near the r point (with q ^ 0) is important for electronic 
energy spectrum. 



Finally, we give two examples of lattice deformation along a line in graphene as 
shown in Fig.[5ja) and (b) whose field B% (r) can polarize the pseudospin. In Fig.^a), 
we modify the hopping integral at the dotted line of y = 0, so that #7o(r) ^ at 
y = and <ho( r ) = ^7o( r ) = °- From Eq. (pTT) . we obtain v F A q (y) = (<57o(y),0). 
In Fig. E{b), the hopping integral along the dotted line at x = (precisely, x = 0+ 
or x = 0_) is changed as #7o(r) = and #7o(r) = #7o(r) 7^ at x = 0. We 
obtain vfA (1 (x) = (— S^q (x), 0) for this case. Figure [5^a) and (b) correspond to 
the generation of the zigzag edge and the armchair edge in the limit of a strong 
perturbation, respectively, if we take = 70 and 5jq = = 70 for the hopping 
integral at the two lines. The behavior of the electronic structure for the two cases are 
completely different from each other, which we show by calculating the deformation- 
induced magnetic field. In the case of Fig. [5fa), we have a finite deformation-induced 
magnetic field Bz(r) 7^ near the line at y = 0. The deformation-induced magnetic 
field Bz(y) is negative at y < as illustrated by (g> and is positive at y > as 
0. On the other hand, in the case of Fig. EJb), the deformation- induced magnetic 



V • u LO (r) / 0, V x u LO (r) = 0. 



(74) 



5 q (r)/0, V-A q (r) = 0. 



(75) 



§5. Edge States of graphene 
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field is zero. The gauge field which gives zero magnetic field can be removed from 
Hamiltonian by a gauge transformation, which is discussed in the previous section. 
The deformation-induced magnetic field accounts for the presence of so-called edge 
states only at the zigzag edge, which is shown as below. As we show in § [IJ the edge 
state is a localized wavefunction near the zigzag edge in which the pseudospin of the 
edge state is perfectly polarized. This means that the amplitude of the wavefunction 
has a value only on either A or B atoms. 

(a) B z {r)^0 (b) B z {r) = 




Fig. 5. Local deformation along a line and corresponding deformation-induced gauge field A q (r) 
are shown by arrows in (a) zigzag and (b) armchair edge. In (a), A q (r) gives a finite deformation- 
induced magnetic field -B q (r) (flux) illustrated by and ® while no deformation- induced mag- 
netic field is present in (b). 

We will derive the edge states from the effective-mass model with the gauge field 
of Fig. Efct). 10 ) Before going into the details, let us first outline the story. We will 
show that there are localized pseudospin-polarized states in the energy spectrum and 
the energy dispersion appears at p x > shown as the two solid lines in Fig. E^b). 
The velocity for the energy dispersion becomes small with increasing the gauge field 
and it becomes zero when the gauge field is sufficiently strong which corresponds 
to the zigzag edge. This result of the effective-mass model can reproduce the result 
of the tight-binding (TB) lattice model as shown in Fig. E|c),(d) and (e), 10 ) where 
the TB lattice model is defined by an adiabatic parameter c as ^7o(y = 0) = C70 in 
Eq. ([27]) . Here c = and c = 1 correspond to no deformation and the zigzag edge, 
respectively. 

We assume that A q (r) of Fig. [5{a) is quite localized within \y\ < £ g , that is, 
Ai(r) = (A§(y),0) and A%(y) = for \y\ > £ g in Eq. flU, where £ g is a length of 
the order of lattice spacing. We parameterize the localized energy eigenstate as 

^|'(r) = N' eMik x x)e- G ^ , (76) 

where N' is a normalization constant. The pseudospin polarization is represented by 
g(y), and the wave vector k x is a good quantum number because of the translational 
symmetry along the x-axis. In the direction of y, we assume a localized nature of 
the wavefunction which was obtained by TB calculation. Putting Eq. (|76p to the 
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Fig. 6. (a) Band structure (k x ,E) around the K' point of the effective- mass model without the 
deformation-induced gauge field and (b) that with the gauge field along the a;-axis. The shaded 
regions in (a) and (b) represent the spectrum of extended states for different k y . The (two- 
dimensional) dispersion relation around the K' point for graphene without the deformation 
gives two cones (representing the linear k dispersion) whose apex is the K' point as is shown in 
the inset of (a). The band structure of (a) is a one-dimensional projection (onto k y plane) of 
the linear k dispersion. In (b), the energy dispersion for the localized state is represented by 
solid lines. When the bonds at y — become disconnected (c = 1) and zigzag edges appear, 
the energy eigenvalues of the localized state converge to zero (E — > 0) for any k x (> 0) and 
forms a flat energy band of edge state. Using the TB model, we plot the band structure for 
an undeformed graphene in (c), for a zigzag ribbon in (e), and for a graphene with weakened 
hopping integral for the C-C bonds at y = in (d). k x a in (c), (d) and (e) is the wave vector 
parallel to the x-axis and k x a — 2ir/3 is the K' point. 



energy eigenequation, Ti. K ipg (r) = Eip^ (r), we obtain 

p x - A%{y) - hj- (G(y)+g(y)) = -^e +2 ^\ 

dy v F (7?) 

Px - A*(y) + h— (G(y) - g(y))) = e~ 2 ^ . 

dy v F 



By summing and subtracting the both sides of Eq. (|77p . we rewrite the energy 
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eigenequation as 

p x - Al{y) - = -— cosh(2 5 (y)), 

dy v F 

dG(y) E , / w 

ay vf 

Since we are considering a localized solution around the x-axis, we assume G(y) ~ 
\y\/£, where £ (> 0) is the localization length. When G(y) ~ |y|/£, the solution of 
the second equation of Eq. ([78]) is given by 

sW H M (»>«. (79) 

The functions and g(y) are schematically shown in Fig. [FJ^a) and (b), respec- 
tively. The sign of g(y) changes across the x-axis; the sign change of g(y) means that 
the pseudospin direction changes at the x-axis. The change of pseudospin is induced 
by the gauge field A%(y). To see this, we integrate the first equation of Eq. (|78p from 
y = —^ g to £ g , and we get the following relation 

-°f±dy = - A%{y)dy. (80) 

We have neglected other terms, since they are proportional to £ g and become zero 
in the limit of £ g = 0. By putting Eq. ([79]) to Eq. ([80]), we find 

- sinh ^(^j = ^j^A^ y )dy. (81) 

When the right-hand side of Eq. (|8ip is large, we obtain from Eq. (|79p that g(y) 3> 
for y < —£g and g(y) <C for y > £ g . In this case, the localized state is a pseudospin- 
up state iPe(t) oc *(1, 0) for y < — £ 5 and a pseudospin-down state iP%(t) oc *(0, 1) for 
y > Cg- Thus, a strong gauge field at the x-axis makes pseudospin-polarized localized 
states. Since the polarization of the pseudospin means that the wave function has 
amplitude only A (or B) atoms, this result is consistent with the result by the TB 
model for the edge state. 9 ) We understand that the edge state is a pseudospin 
polarized and localized state in the real space which is induced by the deformation- 
induced magnetic field of the A q (r) along the x-axis. 

Let us now calculate E and £. To this end, we use the first equation of Eq. ([75]) 
for \y\ > £ g and obtain 

cosh - / A*(y)dy 




Moreover, using Eq. ([HT]) . we get 

h . ( i f <9 




A%{y)dy . (83) 
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Fig. 7. (a) The amplitude of the wave function, exp(— G(y)), of a localized state whose localization 
length is £. (b) The pseudo-spin modulation part, g(y). From Eq. (|79[) , g(y) is a constant for 
\y\ > £s> an< l abruptly changes across the x-axis (y — 0). 

In addition to this localized state, there is another localized state for the same p x 
with the same ^ but with an opposite sign of E. This results from a particle-hole 
symmetry of the Hamiltonian; a z 7i K a z = —7i K . By the particle-hole symmetry 
operation, the wave function is transformed as iP_ e (y) = o z ^, (r). 

The normalization condition of the wave function requires that £ should be 
positive, which restricts the value of p x to p x > in Eq. (|83p . In fact, when A x (y) is 
positive, Eq. ([83j) means that the localized states appear only at p x > around the 
K' point. This is the reason why the localized states appear in the energy spectrum 
only in one side around the K' point as shown in Fig.[6^b). On the other hand, in the 
case of the K point the Hamiltonian is expressed by Eq. (|32|) . Because of the different 
signs in front of A q (r) of Eqs. (|32p and (|33p . a similar argument gives p x < for the 
edge states around the K point. Furthermore, when (1/ft) A x (y)dy 3> 0, E in 
Eq. (|82p becomes zero. The zero energy eigenvalue between the K and K' points in the 
band structure corresponds to the flat energy band of the edge state. 9 ) When A x (y) 
is negative (A x (y) <C 0), a flat energy band appears in the opposite side: p x < 
around the K' point and p x > around the K point. This condition, A x (y) -C 0, 
corresponds to the Klein's edges 31 ) which are obtained by removing A or B atoms 
out of the zigzag edges having A or B atoms. Calculated result by the TB model 
with the Klein's edges 31 ) is thus analytically explained by A q (r). 

There are also extended states in addition to the edge states. The energy dis- 
persion relation of extended states can be obtained as (E/vf) 2 = P x +Py, by setting 
G(y) ~ —ip y y/h (\y\ > £ g ) in Eqs. ([76|) and ([78]) where p y is a real number. The 
calculated energy bands are given by \E\ > vf\p x \, shown as a shaded region in 
Fig. [6]^b). It also agrees well with the TB calculation shown in Fig. [6^e). Since 
1 — tanh 2 x = 1/ cosh 2 x, we see that the energy dispersion relation of the localized 
state between E and £ becomes 



which is the same as the linear dispersion relation (E/vf) 2 = P x +P 2 if one replaces 
p y with i(h/^). Thus One can then regard the localized state as a state with a 
complex wavenumber. 




(84) 
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We have shown that three basic properties of the edge states, i.e., the pseudospin 
polarization, the dependence on the momentum, and the flat energy band, obtained 
previously by the TB model, 9 ) can be explained analytically in terms of the gauge 
field. In order to compare the present theory with the TB model quantitatively, we 
have performed a TB calculation for the geometry of Fig. 0a) with setting = 070- 
In Figs. E^c), Eld) and[6^e), we plot the band structure for c = 0, c = 1/2, and c = 1. 
Comparing these figures with the results of the continuous model, one can find an 
exact correspondence between the TB model and continuous model. Moreover, we 
analytically find 10 ) 

v F cosh(-ln(l - c)) ^ 
§ = Px tanh(- ln(l - c)) + 0(lp 2 x /h), 

for localized states around the K' point. Thus, by comparing Eq. (|85j) with Eqs. (1821) 
and (|83p . we conclude that the TB model and the continuous model agree with each 
other near the K' point (p x l/fr <C 1), by the following relationship 

1 & 

A x (y)dy = -ln(l-c). (86) 



n 



The right-hand side diverges when c — ► 1, which reproduces the flat energy band 
(E -> 0) in Eq. (|S2J) and gives £/h = p' 1 in Eq. ([83]). When c -> 0, we have 
^/_| Ax(y)dy — c, which confirms Eq. (|3ip that is derived by assuming a weak 

perturbation for #7o(r).E3 

By considering the edge state using the effective-mass model, we found that the 
deformation- induced gauge field (A q (r)) and magnetic field (5 q (r)) explain basic 
properties of the edge state. It is summarized as follows: 

(1) The gauge field can generate the edge state in energy spectrum, depending on the 
gauge field direction. Let ey the unit vector along the edge and ^(r) = A q (r) • ey, 
the edge state appears if the gauge field has a component parallel to the edge: 
-<4f|( r ) ^ 0. This explains the presence (absence) of the edge state near the zigzag 
(armchair) edge. 

(2) The edge states are pseudospin polarized states. The direction of the pseudospin 
polarization is determined by the direction or sign of I? q (r). 

(3) The direction of the gauge field, namely, ^4 q (r) > or ^4 q (r) < 0, is vital for the 
edge states, since it determines the energy dispersion and wave vectors which allow 
the edge states. 

(4) The flat energy band of edge states at zero energy is obtained by the limit; 
A^(r) -> ±00. 



*' It is interesting to note that -B?(r) for the edge state corresponds to an enormous magnetic 
field ~ f0 5 T at the zigzag edge. 10 -*' 32 - 1 Thus a uniform external magnetic field has little effect on 
the edge state, compared with the deformation. 
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§6. Discussion and Summary 

By formulating the effective-mass equation for a graphene with a lattice defor- 
mation, we have shown that the deformation can be represented by the deformation- 
induced gauge field, A q (r). We formulate the el-ph interactions and represents the 
shape of edge by A q (r). The appearance of the gauge field in the effective- mass 
equation is reasonable because the Feynman diagram for the el-ph interaction is 
basically the same as that for the electromagnetic interaction of A(r). The only dif- 
ference between the fields A q (r) and A(r) is that A q (r) does not break time-reversal 
symmetry as a whole system while A q (r) breaks time-reversal symmetry locally in 
k-space. Thus, we think that a time-reversal symmetric gauge field is useful not only 
for graphene system but also other lattice systems that have an internal degree of 
freedom like the pseudospin. 

Here we would like to mention the extension of the effective-mass model for 
the edge states. It is known that the Coulomb interaction makes the real spin 
(not pseudospin) of the edge states polarized. 9 ** We considered the effective-mass 
model for the Hubbard model with a mean field approximation, and found that the 
Hubbard term appears as a mass term whose sign depends on the spin of the edge 
states. The mass term shifts the energy of the edge states up or down relative to 
E = depending on the sign of the mass so that the ground state exhibits the spin 
polarized state. The gauge field and the mass term of the effective-mass model give 
rise to a parity anomaly phenomena. We have shown that ferromagnetism of the 
edge states closely related to the parity anomaly. 32 ) 

In summary, the lattice deformation of graphene is modeled in the effective- 
mass theory by the deformation-induced gauge field which keeps the time-reversal 
symmetry but breaks the special symmetry defined within each valley. The formalism 
of the gauge field and the gauge symmetry is useful for understanding anomalous 
physical properties of graphene. 
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